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Abstract 

We have studied the electronic properties of the ladder compound a'-NaV205 
adopting a joint experimental and theoretical approach. The momentum- 
dependent loss function was measured using electron energy-loss spectroscopy 
in transmission. The optical conductivity derived from the loss function by 
a Kramers-Kronig analysis agrees well with our results from LSDA+C/ band- 
structure calculations upon application of an antiferromagnetic alignment of 
the V 'idxy spins along the legs and an on-site Coulomb interaction U of 
between 2 and 3 eV. The decomposition of the calculated optical conduc- 
tivity into contributions from transitions between selected energy regions of 
the DOS reveals the origin of the observed anisotropy of the optical conduc- 
tivity. In addition, we have investigated the plasmon excitations related to 
transitions between the vanadium states within an effective 16 site vanadium 
cluster model. Good agreement between the theoretical and experimental loss 
function was obtained using the hopping parameters derived from the tight 
binding fit to the band-structure and moderate Coulomb interactions between 
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the electrons within the ab plane. 
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I. INTRODUCTION 



In 1996 M. Isobe and Y. Ued J published magnetic susceptibility measurements of a'- 
NaV205 powder samples, proposing the existence of linear antiferromagnetic spin-1/2 chains 
and a possible spin-Peierls transition at a critical temperature, = 34 K. The original pic- 
ture of the charge ordering, in this mixed valence ladder compound was that of alternating 
legs of V^"*" and V^"*" ionsB However, a recent determination of the crystal structure using 
single crystal x-ray diffraction at room temperatureiS yielded only one symmetrically in- 
equivalent vanadium (V^'^^) position. Due to the observed one-dimensional character of the 
spin and charge system at room temperature - which was by that time also confirmed by 
inelastic neutron scattering!'! and angle-resolved photoemission experiments! - the picture 
of linear chains of V-O-V rungs containing a single ci-electron in a molecular orbital-like 
state with antiferromagnetic alignment along the ladders has increasingly won recognition. 
Below Tc charge ordering has been observed in V NMR studies,i but the ordering pattern 
has not yet been defined. Calculations of the electronic ground stateB0 have predicted a 
zig-zag charge orderingalong the ladders and only lately, with the support of the most re- 
cent experimental datalHTtJ does the evidence for zig-zag order appear to outweigh that for 
more complicated order patterns.^ Nevertheless, there remain a number of important, open 
questions such as the nature of the correlations between the charges and spins which are the 
driving force for the observed charge ordering and which make the system insulating. 

The electronic structure of a'-NaV205 has been intensively studied experimentally by 
means of optical absorptionlll and reflect ivity,llllli as well as theoretically within the frame- 
work of bandstructure calculationg^EjEll and exact diagonalization techniques.cS^E^ Never- 
theless, there is a controversial discussion about the origin of the electronic excitations in 
this system. For example, concerning the absorption peak occurring at about l^V in optics, 
transitions between bonding and antibonding combinations of V Sd^y states,i'0ii transitions 
between vanadium states of different svmmetryil as well as on-site d — d transitions between 
crystal-field-split vanadium 3d statesli^ have been proposed as the origin of this feature. 
Furthermore, the information from the experiments still appears to be insufficient to enable 
the definition of a unique parameter set for the description of the electronic structure of 
a'-NaV205. In particular, independent of which theoretical model was used, significantly 
different values for the on-site Coulomb repulsion U, the inter-site Coulomb repulsion Vxy 
between electrons on legs of adjacent ladders and the hopping parameter t^y have been dis- 
cussed in the literature (for a comparison see Table |). Since the charge ordering as well as 
the charge transport properties strongly depend on the energies described by these parame- 
ters, a more exact determination of their values could lead to a better understanding of the 
electronic properties of a'-NaV205. 

In this contribution, we present a joint experimental and theoretical investigation of the 
electronic excitations and their momentum dependence in a'-NaV205 measured using high 
resolution electron energy-loss spectroscopy in transmission. The comparison of the data 
with models based both on the LSDA+U formalism, as well as a cluster approach, enable the 
construction of a consistent theoretical description of the electronic structure of a'-NaV205 
and result in the determination of a double-checked parameter set. The paper is organized 
as follows. In Sec. II the experimental and theoretical methods are introduced. Section III 
contains the presentation and discussion of the experimental and theoretical results and is 
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split into sub-sections dealing with the different aspects of the data. Lastly, Sec. IV is a 
summary. 

II. METHODOLOGY 
A. Experiment 

1. Samples 

Single crystals of a'-NaV205 were grown from the melt. The detailed procedure is 
described elsewhere.0 a'-NaV205 crystallizes in an orthorhombic unit cell with lattice con- 
stants a = 11.318 A, b = 3.611 A and c = 4.797 A.i The vanadium atoms and their five 
nearest neighbor oxygen atoms build slightly distorted VO5 pyramids which are connected 
by their corners along the b direction as well as by their corners and edges along the a 
direction, as shown in Fig. |I|a. The resulting pyramid layers in the ab-plane are separated 
by sodium atoms. The weak interactions between the vanadium oxide layers and the sodium 
layers are the reason for the good cleavage behavior of this system along the plane perpen- 
dicular to the c direction. The projection of the vanadium and the oxygen atoms of a single 
pyramid layer onto the ab plane, excluding the oxygen atoms on the apices of the pyramids, 
delivers the atomic configuration as illustrated in Fig. |I|b. This approach emphasizes the 
ladder structure in which adjacent ladders are shifted with respect to each other by half of 
the lattice parameter b in the b direction (see dashed lines in Fig. |l]b). 

For the measurements using electron energy-loss spectroscopy in transmission, thin films 
of about 1000 A thickness were cut from the single crystals with a diamond knife using 
an ultramicrotome. Because of the good cleavage behavior, the crystallinity remains con- 
served after cutting parallel to the ab plane. The high quality and orientation of the single 
crystalline samples were checked by in situ electron diffraction. 

2. EELS in transmission 

EELS in transmission with a primary beam energy of 170 keV was performed on free 
standing films at room temperature (for experimental details see Ref. ^). The energy and 
momentum transfer (q) resolution were chosen to be 110 meV and 0.05 A~^ for q < 0.4 A~^, 
and 160 meV and 0.06 A~^ for q > 0.4 A~^, in order to compensate for the decrease of the 
cross section at higher momentum transfer. 

EELS in transmission provides us with the momentum and energy dependent loss func- 
tion lm(- l/£:(q, w)), from which, by means of the Kramers- Kronig relations, the real part 
of the negative inverse dielectric function and thus all the optical properties such as, for 
example, the optical conductivity cr{uj) can be calculated. For small momentum transfer 
only dipole transitions are allowed, and for the limit q = the transition matrix elements 
are the same as in optics. For the Kramers- Kronig analysis the loss spectra closest to the 
optical limit (q = 0.1 A~^) were used in order to derive the optical conductivity spectra. 
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B. Theory 



1. Band structure calculations 



The band structure of a'-NaV205 was calculated self-consistently using the scalar- 
relativistic LMTO methodS in the atomic sphere approximation with the combined cor- 
rection taken into account (ASA+CC). The von Barth-Hedin parameterizationil was used 
for the exchange-correlation potential constructed in the local spin density approximation. 
The Brillouin zone (BZ) integrations in the self-consistency loop were performed using the 
improved tetrahedron method.0 

The experimental lattice parameters and atomic positions as determined for the high 
temperature phase (Ref. ^) were used in the calculations. The angular momentum expansion 
of the basis functions included / = 3 for vanadium and sodium and I = 2 for oxygen and the 
empty spheres.^ The O 3d and the V 4/ states were included in the basis, as they give a 
significant contribution to the dipole matrix elements of the optical transitions. 

In order to account for the strong electronic correlation at the V sites, we used the 
LSDA+[/ method in the band structure calculations!!, which has been shown to be very 
helpful for the description of the electronic structure of transition metal oxides, in which the 
3d orbitals hybridize quite strongly with the oxygen 2p orbitals (for a review see Ref. pB] ). 
In this method a Hubbard-like term is added to the LSDA total energy functional: 

^LDA+C/ ^ ^LSDA _ ^dc ^ ^-^^ 

where E^^^^ is the LSDA energy functional, E'^ takes into account the on-site Coulomb and 
exchange interactions, and E'^'^ is necessary to avoid the double counting of the averaged 
Coulomb and exchange interactions already included in iv;LSDA non-spherical contribu- 
tions to the on-site Coulomb and exchange integrals U and J are neglected, the E^ term 
can be written as: 

= {U - J)ni^nj^ + ^Y^Uni„nj_„, (2) 

where is the occupancy of the i-th localized orbital with the spin projection a. The 
localized orbitals used in are constructed in such a way that they diagonalize the charge 
density matrix ra^-. Then, the expression for the orbital dependent one-electron potential 
corresponding to (|l]) is given by: 

yL^A^u ^ ^LSDA ^ j2iU - J)(i - na^a >< ^a\ (3) 

where \ia >< ia\ is the projector onto the localized orbital. 

The absorptive part of the optical conductivity tensor components was computed from 
the LMTO energy bands and eigenvectors on the basis of the linear-response expressions,@H2l 
whereas the dispersive part was obtained via the Kramers-Kronig transformation. The finite- 
lifetime effects and the experimental resolution were simulated by broadening the calculated 
dielectric tensor spectra with a Lorentzian width of 0.2 eV. 

As mentioned in the introduction, a'-NaV205 can be considered as composed of linear 
chains of V-O-V rungs containing a single d- electron in a molecular orbital-like state with 
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antiferromagnetic alignment along the ladders. For the LSDA+f/ calculations presented 
here, initially an antiferromagnetic (AFM) order of each of the V spins along the leg of the 
ladder is assumed. We have also carried out a brief analysis of the influence of deviations from 
this AFM order upon the calculated DOS and optical conductivities by means of spin-waves. 
The band structure calculations for non-coUinear magnetic structures including constrained 



moment directions were performed using the formalism developed in the Refs. ^ and |39 . 
The magnetization direction in the atomic sphere surrounding the site t + Rj, where t is its 
position in the unit cell and Rj is a real space lattice vector, is defined by two polar angles 
9ti and (pu- If Ou does not depend on Rj and (pu is given by 

(j)u = (f)t + u- Ri, (4) 

{(pt determines the relative orientation of the magnetization directions for different sublat- 
tices), the calculations can be performed without an enlargement of the unit cell for an 
arbitrary vector u, which defines a spiral magnetic structure in real space. (In the present 
paper we adopted the notation u instead of the usually used character q to distinguish the 
vector of a spiral structure from the momentum transfer vector in the EELS measurements_L. 
The details of the U calculations for spiral magnetic structures will be given elsewhere.EJ 
Here we only mention that in this case the effective orbital dependent potential is still given 
by equation (^) but with the spin part of a localized orbital \ia > centred at a given atomic 
site being defined in the local coordinate system with the z axis parallel to the magnetization 
direction in the corresponding atomic sphere. 



2. Cluster calculations 



We have studied the dynamic dielectric response using a quarter-filled t-J-V model 
H = - h (4,<xCi,a + H.c.) + (Si ■ - i 



4 " ^ 



jTliTlj 



(5) 



where c- = c- ^.(l — rii-a) are constrained electron creation operators, ni = Y.a (^i,a'^i,a- is the 
occupation-number operator, and Si denotes the spin-i operator at site i. The expression 
{i,j) denotes the summation over all pairs of nearest neighbors. The hopping parameters 
tij, intersite Coulomb interactions Vij, and exchange interactions Jij are defined in Fig. ^ 

The loss function measured in EELS experiments is directly proportional to the dynamic 
density-density correlation function.Ell By including the long-range Coulomb interaction in 
the model within a random-phase approximation (RPA) one finds for the loss function 



L{uj, q) = Im 



1 + ^qXp(^>q). 



where 



Auit 



(6) 



(7) 
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is the response function at zero temperature for the short-range interaction model described 
by Eq. (^). x° depends on the energy loss uj and momentum transfer q. |0) is the ground 
state, pq denotes the Fourier transform of n^, and fq = N / {eQerVc^) is the long-range 
Coulomb potential with unit cell volume v. Furthermore, N is the number of electrons per 
unit cell, eo is the permittivity, and er is the real part of the dielectric function. 

We evaluated Eq. by direct diagonalization using the standard Lanczos algorithm.Ei 
As this method is limited to small clusters, the maximum cluster size in our calculations was 
restricted to 16 vanadium sites. However, calculations for a cluster consisting of two adjacent 
ladders, each containing four rungs, showed strong finite-size effects. We have solved this 
problem by using two different clusters for momentum transfer parallel to a and b direction, 
as denoted in Fig. ^ by the filled circles. In the a direction we chose an array of double 
rungs extending over four adjacent ladders, whereas in the b direction the cluster consisted 
of a single ladder with eight rungs. The separation into these two clusters is justified by 
the small value of the interladder hopping amplitude t^y = 0.012 eV (details will be given 
elsewhereS) . This value of t^:^ has also been found by previous LDA calculation^, as well as 
from a tight binding fit to our own LDA results.0 Additional evidence for a small txy follows 
from the weak magnetic dispersion along the a direction as observed in neutron scattering 
experiments.!^ 

For the calculation of the loss function we chose open boundary conditions, and for 
momentum transfer parallel to the a direction we used renormalized intersite Coulomb 
interactions Va = Va + Vb and Vb = 2Vb. These values follow from a straightforward analysis 
of the infiuence of adjacent rungs of the same ladder on the electronic states.^! Furthermore, 
in the case of open boundary conditions one has to make sure that electrons on the edges 
of the cluster are still embedded in the local Coulomb potential that results from a zig-zag 
charge ordered state. Thus, sites on the edge of the cluster for momentum transfer parallel 
to the a direction are assigned an additional on-site energy Vxy. By analogy, sites on the 
edges of the cluster for q parallel to the b direction need an additional on-site energy Vb, if 
they are not occupied in a zig-zag charge ordered state. 



III. RESULTS AND DISCUSSION 

A. Loss function with small momentum transfer 

We have measured the energy dependent loss functions for different momentum transfers 
q parallel to the crystallographic a and b direction in an energy range between 0.5 eV and 
70 eV. Due to the contribution of the elastic line and surface losses it is not possible to 
measure at zero momentum transfer but close to the optical limit (q = 0.1 A~^). In Fig. |^ 
the loss functions with momentum transfer q = 0.1 A~^ parallel to the a (solid line) and 
the b direction (dashed line) are shown. Both spectra are dominated by a broad feature at 
around 23 eV which represents the volume plasmon - a collective excitation of the valence 
electrons. While the feature at 50 eV results from local excitations between the V 3p and 
the V 3d states, the features emerging below the volume plasmon energy can be assigned to 
transitions of the valence electrons into the conduction band. As Fig. ^ shows, an anisotropic 
behavior is only visible in the low energy regime (< 22 eV). 
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In order to compare the spectra with results derived from our LSDA+f/ calculations, 
we have calculated the optical conductivity, u, by performing a Kramers-Kronig analysis 
(KKA). The comparison with the optical conductivity of the structurally related V2O5 (d° 
configuration) obtained with the same method,^ shows that above ~5 eV the spectra of both 
compounds look very similar (not shown). In Ref. the anisotropy of the loss function and 
of the optical conductivity of V2O5 are discussed in detail. As a result, the features above 
5 eV can be assigned to transitions from occupied bands with mainly oxygen character into 
unoccupied bands which are dominated by the V 3d states. From the resemblance of the 
spectra and the related crystal structures of both compounds we conclude that the same 
holds for a'-NaV205. Thus, in this work we focus on the energy range below 5 eV where 
the loss functions and, consequently, the optical conductivities of a'-NaV205 differ strongly 
from those of V2O5. 



B. DOS 

In comparison to V2O5, which has a similar network of linked VO5 pyramids, the Na 
atom in a'-NaV205 donates an additional electron to the valence band. Since in V2O5 
the bands formed by the O 2p states are essentially completely filled and the energetically 
lowest lying unoccupied bands are formed by the V Sd^y states,^ the additional electron 
in a'-NaV205 goes into the rather localized V 3dxy states. These states determine, to a 
great extent, the fascinating properties of this compound. Due to the localized nature 
of the V 3d electrons, the Coulomb correlations between them are rather strong and, as 
a consequence, conventional L(S)DA calculations fail to describe properly the electronic 
structure of a'-NaV205. In the unit cell corresponding to the centrosymmetric Pmmn 
structure of the high temperature phase, a non-magnetic metallic solution is obtained, in 
contradiction to the experimental data. One way to overcome this discrepancy is to introduce 
an antiferromagnetic (AFM) order of V magnetic moments along the b direction. This leads 
to a doubling of the unit cell and the opening of an energy gap in the LSDA band structure. 
However, the value of the energy gap is still strongly underestimated. It should, however, 
be borne in mind there is no experimental evidence for the existence of AFM order in 
a'-NaV205 even at low temperatures. 

An efficient way to take the strong correlations in the V 3d shell into account is to use 
the LSDA+f/ approach. However, this method is usually applied to ordered compounds and 
thus an additional assumption has to be made in order to study the electronic structure and 
the optical properties of the high temperature phase of a'-NaV205. Since, according to the 
experimental data, there is no charge ordering above Tc and all V sites are equivalent,! we 
have assumed that the effective one electron potential at each V site is the same and equal 
to the average of the LSDA+[/ potentials for and ions. In other words, the orbital 
dependent potential is calculated using the occupation numbers averaged formally over the 
V'*^ and V^"*" orbital occupation numbers. On application of U, the partially occupied V 
dxy orbitals are the most affected. The average occupation of the majority-spin dxy orbitals 
is close to 0.5 and their energy position remains unchanged as compared to LDA, whereas 
the unoccupied minority-spin dxy orbitals are shifted upwards by f//2. As a result, an 
insulating behavior - with the magnetic moment of 0.5 /i^ per V atom - or 1 /i^ per rung - is 
immediately obtained, which is independent of the kind of magnetic order which is assumed 
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along the b axis. One can say that the on-site Coulomb repulsion effectively suppresses the 
occupation of the d^y orbitals with the opposite spin at each rung. 

The effect of the absence, in reality, of long range AFM order in the high temperature 
phase of a'-NaV205 on the electronic structure and optical properties was modelled by 
performing calculations for spin waves defined by a vector u = (0, Uy, 0) with Uy varying 
from 0.5 (corresponding to AFM order of V magnetic moments along the b direction), to 
at which all V magnetic moments are ferromagnetically aligned. In these calculations all dt 
were set to zero, whereas the (pt were chosen in such a way that the V ions situated at the 
same rung had the same magnetization direction. We found that the relative orientation of 
the magnetization at V sites of neighbouring ladders, which is also determined by (pt, has 
only minor effect on the calculated band structure. 

The calculations were iterated to self-consistency for Uy = 0.5 and then one iteration 
was performed for spin waves with other values of Uy. The comparison of the band energies 
confirmed that the ground state of the system is antiferromagnetic, in accordance with 
previous estimates for the exchange coupling constants between V spins in the b direction. 

The densities of the V 3d and O 2p states, calculated for Ucs = 3 eV and an AFM 
order of the V magnetic moments along the b direction are shown in Fig. ^. In Fig. ^b, the 
density of V 3d states calculated for Uy = 0.3 (dashed lines) is also shown. Here, and in the 
following, we assume that the magnetic moments of the V ions situated at the same rune 
are ferromagnetically ordered. The overall structure of the DOS is similar to that of V2O5.E3 
The states in the energy range from -6.5 to -2 eV are formed mainly by the O 2p states 
(Fig. ||a) with a bonding hybridization with the V 3d states. The V dxz,yz states dominate 
in the 1.8-3 eV range, while the states with V d3z^_r2 and V d^i-y^ character are shifted to 
higher energies due to the comparatively strong hybridization with O 2p states and form the 
upper part (> 3 eV) of the conduction band, shown in Fig. The main differences in the 
structure of the DOS as compared to V2O5 are caused by the additional electrons occupying 
the V dxy bands which are empty in the case of V2O5. Since the change in the occupation of 
the V dxy states is the main source of the differences between the optical spectra and other 
physical properties of these compounds, in the following we will focus our attention on the 
states originating from the V dxy orbitals. 

The narrow peak in the majority-spin V 3d DOS just below the Fermi level (Fig. ^3) 
originates from a combination of the 3dxy orbitals of the two V atoms on the same rung which 
is orthogonal to the 2py orbital of the Or atom in the center of the rung. The corresponding 
antibonding V 3dxy-0^ 2py states are responsible for the DOS peak at 1 eV. In terms of the 
effective V-V hopping these two peaks result from the bonding and antibonding combination 
of the V 3dxy states, respectively, and the energy separation between them is determined 
by ta = 0.38 eV acting across the rung. Here, for the effective V-V hopping terms we use 
the values obtained from a tight binding fit to the LDA band structure calculated with the 
LMTO method^ which were found to be very close to the hopping parameters determined 
in Ref. |^. The minority spin dixy states of bonding character are shifted to higher energy by 
the effective Coulomb repulsion and lie above the antibonding d^xy peak. While the energy 
of the latter is governed by ta and thus does not depend on Ues, the relative position of 
the minority spin dixy states does depend on the exact value of Ues- Finally, the peak at 
~2.5 eV arising from the antibonding dixy states is almost lost in the contributions from the 
other states of the conduction band, mainly formed by the remaining V 3d states hybridized 
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with O 2p states. 

It should be noted that for both spin directions the width of the bonding d^y peaks 
in the LSDA+t/ calculations is much smaller than 0.7 eV, the value obtained from LDA 
calculations, because the relatively strong hybridization along the b direction governed by 
th = 0.17 eV is suppressed due to the AFM order. However, as Uy - which defines the 
magnetic structure along the b direction - decreases, the width of the peaks increases and 
reaches the LDA value when the V magnetic moments are ordered ferromagnetically. As 
for the antibonding dxy states, the dispersion of the corresponding bands along the T-Y 
direction is much weaker already in ldaS and the width of their DOS does not depend on 

Uy. 

C. Optical conductivity 

In order to verify our theoretical results we turn our attention to the optical data. In 
Fig. 1^, the optical conductivity, derived by means of a Kramers-Kronig analysis of the EELS 
spectra measured at low momentum transfer (q = 0.1 A), aligned either parallel to the 
crystallographic a or b direction (Figs. |^a and P^) is shown, together with a^x and ayy 
calculated for Uy = 0.5 using the LSDA {Ucs = 0) and LSDA+[/ with Ucs = 2 and 3 eV 
(Figs. ^ and |^cQ. Our experimental spectra agree well with the previous results derived 
from reflectivityO'lli and absorptionlll measurements of a'-NaV205. Comparing axx with the 
q|| a experimental data (Figs. ^ and ^c), one can see that for this polarization the theoretical 
spectra are only weakly affected by U and thus that the LSDA curve agrees well with the 
experimental one. The calculated position of the low energy peak in centred at 1 eV 
does not depend on U, and coincides almost perfectly with the position of the corresponding 
experimental maximum. However, the calculations overestimate the magnitude of the peak 
and fail to reproduce accurately its asymmetric shape. The peak corresponding to the 
experimental feature at 3.4 eV lies at slightly lower energy in the LSDA calculations and 
shifts upwards to ~3.7 eV with increasing U. 

For (Tyy (Figs. Pd and ^), the situation is quite different. In this case the LSDA calcula- 
tions result in an optical gap of 0.2 eV, which is significantly smaller than the experimental 
value (0.7 eV). In contrast to axx, the position of the low energy peak in ayy is sensitive to 
U, and gets close to the experimental value for Ues = 3 eV. A less strong but still notice- 
able dependence on U can be found in the energy range between 2.5 and 4 eV. Comparing 
the theoretical spectra calculated with different values of f/eg to the experimental ones, one 
can conclude that the best overall agreement between the theory and the experiment is 
achieved for Ues lying in the range between 2 and 3 eV. This range for Ues derived from 
the LSDA+U calculations gives an important guideline for the U value taken in the cluster 
model calculations which will be presented later in the paper. 

The observed correspondence between the calculated and the experimental spectra allows 
us to draw additional information from the detailed analysis of the theoretical spectra. In 
order to understand better the origins of the features of the optical conductivity and the 
reasons for their different dependence on the value of U, we analyzed the contributions to 
the optical conductivity originating from the interband transitions between the initial and 
final states of V 3dxy character. In Fig. ^ we show the results of a decomposition of axx and 
ayy, calculated with Ucs = 3 eV. In each case, the thin solid line marks the sum of all the 
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contributions - i.e. the total optical conductivity. Figs. |]a and ^jc show the calculations for 
Gxxi and Figs. Pd and ||d for Oyy. 

Starting with the top two panels (Figs, ^a and the areas shaded grey indicate spectral 
weight connected with transitions starting from the occupied V "id^y bands. Thus, it is clear 
that the lowest lying peaks in both o^x and Oyy are derived from transitions involving V "idxy 
initial states. The higher lying spectral weight, however, involves transitions starting from 
O 2p states (i.e. the majority of the total o at higher energy is unshaded). 

We turn now to an analysis of the final states involved (see Figs, ^c and ^). In the 
lower panels of the figure we code the various final state characters with different shading: 
(i) the antibonding V "id^xy final states are depicted as light grey; (ii) the bonding V 'id\^xy 
final states are depicted as dark grey and (iii) the hatched areas show the contributions from 
interband transitions into higher lying V Sc? states. 

On comparing the top and bottom panels of Fig. ^ we can therefore conclude the following 
about the transitions giving rise to the optical conductivity. Firstly, the 1 eV peak in Oxx 
is due to a transition between bonding and antibonding hybrids involving the V dxy levels 
of equal spin direction. As was mentioned above, the bonding-antibonding sphtting of dxy 
states is determined only by the hybridization within a rung and is thus not affected by U . 

Secondly, for momentum transfer parallel to the b direction, the matrix elements for 
the bonding-antibonding transitions are zero. In this case, transitions between bonding 
V dxy levels with opposite spin become active. The energy position of these transitions is 
sensitive to the strength of the on-site Coulomb interaction, which naturally explains the 
strong dependence of the position of the low energy peak in Oyy on the value of U observed 
in Fig. 1^. In the language of a many-body approach, the d^xy — ^ d^xy transitions can be 
thought of as excitations which lead to the creation of doubly occupied rungs.S 

Having dealt with the low-lying features of the optical conductivity, we now turn our 
attention to the energy range between 3 and 4 eV. In this context, the lack of grey shaded 
weight in Figs. |^ and ^ illustrates that the initial states of the transitions responsible for 
a in this energy range have O 2p character. Consideration of Figs. ^ and ^d reveals that 
the peak in Oxx at 3.7 eV and the small hump in Oyy at the same energy are dominated 
by transitions into the d^xy bands (dark grey shaded spectral weight). The initial and final 
states of the transitions giving rise to Oxx between 3 and 4 eV are, in fact, of the same 
character as those which determine the shape of the Oxx spectrum in V2O5 just above the 
absorption threshold.Ei The dipole transitions from the O 2j9 valence band states into the 
antibonding d^^y states are completely suppressed for a polarization parallel to a [oxx^-, while 
for Oyy they play an important role in the formation of the shoulder in this energy range 
(shaded light grey in Fig. ^). 

Finally, we mention that transitions from the occupied dxy bands to that part of the 
conduction band originating from the V dxz,yz states give rise to the non-vanishing intensity 
in the calculated spectra in the energy range between the main features (i.e. between 1.5 and 
3 eV). The experimental a's (see Figs. ^ and P?) confirm the accuracy of the theoretically 
predicted optical conductivity in this regard. 

In order to summarize the detailed information won in this section as regards the char- 
acter of the states involved in the transitions giving rise to both Oxx and Oyy in a concise 
form. Table |I| shows a breakdown including the respective initial and final states 'behind' 
each spectral feature. 
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The level to which the detailed breakdown of the calculated optical conductivity has been 
carried out in Fig. ^ is justified by the fact that the theoretical curves (for example the solid 
line) in Fig. || are able not only to reproduce the main features in a^x and ayy quite well, but 
also the weaker, double peak structure in ayy in the 3-4 eV range in our experimental data 
(see Fig. Qd). We note that the spectra derived from optical measurements did not resolve 
this double-peak feature.lill 

Despite the quite good agreement between the calculated and measured optical conduc- 
tivity on the qualitative level, the absolute magnitude and strongly asymmetric shape of 
the low energy peaks in a^x and ayy are not accurately reproduced by the calculations (see 
Fig. We point out two possible sources for these discrepancies: (i) an oversimplified treat- 
ment of the charge fiuctuations in the high temperature phase using an averaged one-electron 
potential and (ii) the assumption that the V magnetic moments are antiferromagnetically 
aligned along the b direction. Dealing firstly with the issue of the charge fiuctuations, we 
note that a more rigorous treatment of the charge disorder on the V sublattice of a'-NaV205 
is beyond the capabilities of the LSDA+?7 approach in its present formulation. We therefore 
return to this question later in the paper in the context of the cluster calculations which 
will be described in detail in the next but one sub-section. As regards the infiuence of de- 
viations from long range AFM order on the optical spectra, we point out that the strength 
of such effects can be estimated by performing calculations for a spin-wave structure along 
the b-axis. 

The optical conductivity spectra calculated for Ucs = 3 eV and for different values of 
the vector describing the spiral magnetic structure in real space, Uy, are shown in Fig. |^. 
Since the change of the magnetic structure affects the bands formed from the V d^y bonding 
states most strongly (see Fig. the spectral features originating from the transitions in 
which these bands are involved demonstrate the strongest dependence on Uy. In particular, 
the change in shape and intensity of the peak at 1 eV in axx (Fig- ^) is the most noticeable 
result of including spiral spin structures along the b-axis. As can be seen from the figure, 
as Uy decreases, the peak becomes less intensive and its maximum shifts to higher energy. 
Since the final states for the corresponding transitions (i.e. the antibonding V d^y states) 
remain very narrow independent of the nature of the magnetic order along the b direction, 
the increase in the width of the peak reflects the change of the density of the occupied V d^y 
states (i.e the initial states). The shape of the low energy peak in ayy (Fig. ^) does not 
change with increasing Uy but its intensity rapidly decreases as Uy decreases, and vanishes 
completely for = at which point the V magnetic moments order ferromagnetically. The 
drop of intensity of the low lying feature in ayy is caused by the decrease of the weight of 
majority-spin states in the final state wave functions. For Uy = the initial and final states 
are formed by pure majority and minority-spin states, respectively, and the corresponding 
transitions are forbidden (since the relativistic effects were not included in the calculations). 
The spectral features in the energy range between 3 and 4 eV are also influenced by Uy, but 
the changes are less pronounced than for the low lying features. 

The comparison of the calculated optical conductivity with inclusion of the spin spiral 
structures (Fig. ^ with the experimental spectra (Figs. ^ and |]b) indicates that an im- 
proved agreement between the theoretical and experimental magnitude of the low energy 
peaks of axx and (J yy CcLIl be achieved if deviations from pure AFM order along the b-axis are 
taken into account. However, even including spin-waves, the LSDA+U results are unable to 
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fully account for the asymmetric shape of the low energy peaks in the experimental optical 
conductivities. 

To try to overcome this weakness, we have also calculated the optical conductivity for 
different relative orientations of the magnetization at neighbouring ladders (not shown). 
The changes of the calculated spectra were found to be very small, which is in keeping with 
the size of the inter-ladder hopping term {txy = 0.012 eV) which couples V ions belonging 
to neighboring ladders. 

To bring this section on the data from band structure calculations to a close, the 
LSDA+t/ results regarding the optical properties of a'-NaV205 can be summarized as fol- 
lows. 

(a) The optical conductivities axx and ayy calculated with Ucs in the range between 2 
and 3 eV are able to reproduce the main features of the experimental spectra, as well as 
the observed optical anisotropy quite well. This fact indicates the relevance and accuracy of 
the parameters determined both from the fit to the band structure (the transfer integrals, 
or t's) or from the comparison with the experiment (U). 

(b) The analysis of the contributions to a coming from interband transitions with dif- 
ferent initial and final states allows us to determine the electronic states responsible for the 
formation of the peaks in the experimental optical conductivity. Table y summarizes this 
information for each spectral feature in a^x and ayy. 

(c) However, the remaining discrepancy between the theory and the experiment as re- 
gards the shape of the low energy peaks (even despite going beyond the AFM approximation 
for the spin order along b) indicates that the subtle details of the electronic structure of the 
high temperature phase of a'-NaV205 are beyond the approximations inherent to the LDA, 
and thus are most likely sensitive to effects more suitably treated within the framework of 
models in which electronic correlation is dealt with at a more fundamental level. 

Before going on to the t — J — V model cluster calculations presented in the last sub- 
section of the results and discussion part of the paper, the next sub-section deals with the 
first EELS spectra of a'-NaV205 recorded beyond the optical limit. Here we exploit the 
strength of EELS in transmission as a measure of the bulk optical properties (in this case 
of a'-NaV205) for finite momenta, thus giving additional insight into the nature of the 
two-particle excitations in this complex system. 

D. Momentum dependent loss functions 

In the left panels of Fig. |^, we show the EELS loss function of Q;'-NaV205 recorded for 
different values of q, either parallel to the crystallographic a (Fig. |^a) or b (Fig. ^) direction. 
The spectra are normalized at higher energy (ca. 7 eV) where the shape of the spectra does 
not change either with the absolute value or with the direction of the momentum transfer 
(see also Fig. For clarity the spectra are incrementally off-set in the y-direction. 

For q parallel to both the a and b directions the spectra show strong intensity between 
the spectral onset at 0.7 eV up to 2 eV, as well as additional weaker features located between 
3 and 4.5 eV. For small q the corresponding plasmon excitations are related to the same 
dipole transitions which are responsible for the peaks at around 1 and 3.5 eV in the optical 
conductivities whose origin was discussed in detail above. In the following, we concentrate 
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our attention on the fine structure of tlie low energy peaks, as well as on the dependence of 
the shape and position of the peaks on the increasing momentum transfer. 

First, we focus on the features below 2 eV. The intensities in the a direction are about 
two times larger than those with the same momentum transfer in the b direction. For q||a, 
the maximum lies at 1.55 eV compared with 1.45 eV for q||b. At high momentum transfers 
the fine structure of the feature around 1.5 eV in the a direction indicates the presence of 
transitions into more than a single final state. One can distinguish between essentially three 
features, whereby the intensity of the component with its maximum at 1.55 eV decreases 
more rapidly with q than the features at 1.2 and 1.7 eV, the latter pair being visible as weak 
shoulders only for high momentum transfers. In the b direction the low energy feature is 
also composed of more than one component - at least at high q. The energetically lowest 
lying feature (at ~ 0.9 eV) - which can be seen as a shoulder at q > 0.4 - looses intensity 
only slowly with increasing q. The second component, which one can relate to the main 
intensity maximum (ca. 1.3 eV) decreases more strongly in intensity at higher q. Lastly, for 
q> 0.5 A~^, a third component is visible at higher energy (ca. 1.7 eV). 

Finally, we turn to the features in the loss function at energies above 3 eV, which, as in 
the data recorded at the low-q limit, are clearly separated from the lower energy features by 
a strong drop of spectral weight between 2 and 3 eV. In the a direction a steep increase of 
spectral weight occurs at 3 eV, which evolves into a maximum at 3.5 eV. With increasing 
momentum transfer the peak position of this feature does not change and the intensity 
decreases slowly. A peak at the same energy position with a similar shape was observed 
for the a as well as for the b direction in the loss functions of V2O5, and is assigned to 
transitions from the O 2p states into the lowest unoccupied V 'idxy states.E3 However, in 
the EELS spectra of a'-NaV205, the shape of the corresponding feature in the b direction 
differs by the presence of a double-peaked structure with maxima at 3.5 eV and 5.1 eV. 
The two peaks merge together with increasing momentum transfer, and for q > 0.7 A^^ 
a further feature becomes visible at ~3.8 eV. Since in electron energy-loss spectroscopy 
the cross section for dipole transitions decreases and that for dipole forbidden transitions 
increases with increasing momentum transfer, the latter 3.8 eV feature is most likely related 
to monopole or quadrupole transitions. 

Having described the experimental loss functions as a function of momentum, we now 
discuss the results of the cluster model calculations aimed at their simulation. 



E. Calculated loss functions 



In Fig. 1^ the results of our cluster calculations (see section [ilB 21 ) are compared with the 



experimental spectra.0 Such cluster calculations require, naturally, input parameters which 
in turn describe the essentials of the physics of the system in question. In this case, for 
the transfer integrals we take the values ta = 0.38 eV, ti, = 0.17 eV, tj:y = 0.012 eV from 
our tight-binding fit to the band structure calculations mentioned earlier. These values are 
identical to those of Ref. |^. A further, vital parameter is the on-site Coulomb repulsion, U. 
Here we take a value of 2.8 eV - a choice which is guided by the range of f/g// = 2 — 3eV 
which came out of the LSDA+f/ calculations presented and discussed above. 

At this stage we can have faith in these parameters for two reasons. Firstly, it is well 
known that LDA does a good job in describing the charge distribution, and hence the hy- 
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bridisation / transfer integrals, even in systems with strong electronic correlation. Secondly, 
the good basic agreement between the optical conductivity calculated within LSDA+t/ with 
experiment (see Sec. |111 C| ) not only lends support to the t's, but underpins the chosen U 
value. 

The exchange interactions were parametrized as Jij = 4:t'^j/U and the values of the 
intersite Coulomb interactions Va = 0.8 eV and Vb = 0.6 eV fulfill the condition for the 
system to be close to a quantum critical point caused by charge ordering, as described 
by Cuoco et al.. il The inter-ladder Coulomb interaction V^y = 0.9 eV has been adjusted 
within the model to obtain the correct peak positions with respect to the experimental loss 
functions. To allow a direct comparison with experiment, the calculated data were broadened 
with a Gaussian function of 0.3 eV width. The unbroadened spectra for q = 0.1 A are also 
shown as vertical lines in Figs. ^3 and As the cluster model contains only V sites, we 
only attempt to simulate the loss function in the low energy range (i.e. for < 3eV). At 
higher energies - as we know from our combined analysis of the experimental and theoretical 
optical conductivities (Sec. [Ill Q ) - excitations from O 2p to V 3d orbitals play a role. Thus 
this spectral weight cannot be described within the present cluster model. 

We find a good agreement between theory and experiment for the low energy part of the 
loss function for q||a, in particular for the data recorded at higher q (see Fig. ^ and ^). 
Specifically, the intensity reduction and increasing width of the experimentally observed 
structure between 1 and 2 eV for q > 0.5 A^^ is well reproduced in the cluster calculation, 
a characteristic which eluded our best efforts within the LSDA+[/ approach. For the calcu- 
lated loss function with q||b (Fig. ^d), the experimental q-dependent reduction in intensity 
is also well reproduced. However, in this case the simulation fails to reproduce accurately the 
large width of the low-energy feature - a shortcoming which becomes increasingly apparent 
for higher q. 

In order to enable a discussion of the origin of the individual features underlying the cal- 
culated spectra, one should note that since a'-NaV205 is close to a charge order transition,ii 
one can describe the nature of the excitations using a zig-zag ordered ground state. In this 
context, we mention that excitations with momentum transfer parallel to the a direction 
lead to a disturbance of the charge ordering as they involve hopping of electrons from one 
side of a rung to the other. Note that this process can also be interpreted as a transition 
from a bonding to an antibonding state of a singly occupied rung, in agreement with the 
results from recent calculations using the Heitler-London modehEl Since the net inter-ladder 
Coulomb interaction, Vxy, remains unchanged after the hopping on the rung in the charge 
ordered state, the excitation energy is dominated by Vb- In addition to the one electron hop- 
ping discussed above, there are also collective processes that involve two or three neighboring 
rungs on different ladders. 

As the unbroadened theoretical results presented in Fig. show, the theoretical spectra 
consist of three main features, each of which can be attributed to a density oscillation 
(plasmon) with energies between 1 and 2 eV. For q||a, the spectral weight shifts from the 
plasmon at 1.6 eV to the one at 1.7 eV with increasing q, which makes it appear as if the 
low lying feature as a whole disperses to higher energy, which is not, in fact, the case. The 
three excitations differ in their degree of delocalization. The excitation at 1.3 eV is related 
to rather delocalized transitions involving more than one electron in the cluster. The higher 
lying plasmons, however, result from more localized processes. 
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For momentum transfer parallel to the b direction, the excitations contributing to the low 
lying spectral weight (see Fig. ^) can be interpreted as transitions from a state in which two 
rungs are singly occupied to one in which one rung is empty and the other doubly occupied. 
Due to the finite hopping amplitude = 0.17 eV, the unoccupied and the occupied rung are 
able to move along the ladder independently from each other. Thus, many more final states 
are available in the b direction than in the a direction. The energies of the transitions into 
these final states depend on the distance between the unoccupied and the doubly occupied 
rung. However, since we can observe only localized excitations by the direct diagonalization 
of small clusters, we are unable to capture fully the dynamics of these excitations and thus 
we miss some spectral weight (in the form of spectral breadth) in the theoretical spectra 
(Fig. ||d) compared to the experimental data (Fig. ||c). This explanation is supported by 
a further reduction of spectral weight in the cluster calculation when we adopt a four-rung 
ladder, in place of the eight rung cluster from which the results are presented here. 

Thus, in contrast to the good description of the strongly localized excitations in the a 
direction described earlier, for q||b we find that the theory significantly underestimates the 
spectral weight, due to finite size effects in the cluster calculation. That these differences 
also persist in the high q regime illustrates that the cluster size imposes limitations even in 
the case where the excitations have shorter wavelength. Nevertheless, it should be pointed 
out that the cluster calculations do provide a qualitative description of both position and 
the q-dependent decrease of intensity of the experimental features. 

The optical conductivity can also be calculated using the same, cluster-based method. 
With the same model parameters we obtain a qualitative agreement with the optical con- 
ductivity, either as derived from EELS or optical measurements^ (not shown). However - 
as was the case for the LSDA-|-[/ data - the asymmetric line shape in a could not be fully 
reproduced in the cluster calculation. We ascribe this to the delocahzed character of some of 
the final states, which at q = (i.e. at infinite wavelength) therefore elude our calculations 
based upon a finite cluster. We note here that a similar type of calculation has been carried 
out recentlyJH and a good agreement between theory and the (q = 0) optical conductivity 
was found, albeit with a rather high value of the inter-ladder hopping integral t^y = 0.15 eV. 
The momentum-dependent EELS data presented here allow a further test of the theoretical 
parameters beyond the optical limit. We have found that adopting t^y = 0.15 eV in our 
t-J-V model calculations yields a significantly worse agreement with the experimental data 
than does the value of t^y = 0.012 eV, which also comes out of the tight binding fit to the 
band structure. Consequently, it appears as if a large value of the interladder hopping is 
not a requirement for an accurate description of the charge excitations in Q;'-NaV205. 

To summarize, one can say that while the small inter-ladder hopping gives a rationale 
for the quasi one-dimensional character of the dynamics of the V 3d electrons; the fact that 
all three inter-site Coulomb interactions {Va, Vb and Vxy) are of similar size is consistent 
with the two-dimensional nature of the charge ordering occurring at low temperatures in 
this system. 

IV. SUMMARY 

In conclusion, we have presented a joint experimental and theoretical investigation of the 
optical properties and collective excitations of a;'-NaV205. From measurements by means 
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of high resolution EELS in transmission, we have derived axx and ayy in the optical limit 
and have also presented the momentum dependence of the loss function with momentum 
transfers parallel to the crystallographic a and b directions. The densities of states and 
optical conductivities were calculated within an LSDA+f/ framework, with the spin order at 
the V sites either described in terms of an antiferromagnetic or spin- wave-like arrangement, 
whereas the momentum dependent loss functions were simulated using a quarter-filled t-J-V 
model based upon a cluster of 16 vanadium sites. 

For q parallel to both the a and b directions, the low energy features of the spectra 
are essentially dispersionless - which is in keeping with the small band width of the unoc- 
cupied states directly above the chemical potential predicted in the LSDA+f/ calculations. 
The comparison between the experimental cr's (measured in the optical limit) and those de- 
rived from the LSDA+f/ calculations yielded a good agreement upon adoption of an on-site 
Coulomb interaction f/ = 2 — 3 eV and an antiferromagnetic ordering of the spins at the V 
sites in the ladders. However, the LDA-based calculations were not able to reproduce the 
width of the experimentally observed structures. 

Based upon the LSDA+f/ results, the differences between a^x and ayy could be analyzed 
in detail and each feature could be characterized according to the nature of the initial and 
final states involved in the underlying optical transitions. For example, it could be shown 
that the optical conductivity between 0.5 and 1.7 eV is dominated by transitions from the 
highest occupied band (bonding combination of V ?>dxy states) into the lowest unoccupied 
band (antibonding V ?idxy-OR 2py states). While for axx such an excitation occurs on a 
single rung, in ayy the same transition can be considered as a hopping process to one of 
the adjacent rungs. The higher lying features in the optical conductivity (above 3 eV) are 
shown to be related to transitions originating in the O 2p valence band manifold. 

Taking the transfer integrals and U value from our LSDA+f/ data (the accuracy of 
which is checked by the comparison with the experimental c's), the t-J-V cluster model was 
successful in describing the q-dependence of the intensity and, in particular the observed 
width of the lowest energy feature at high q for q||a. As the high q data probe shorter 
wavelengths, it is clear that the cluster extension in the a direction is sufficient to describe 
these localized excitations. The situation is different for the data with q||b. Here the 
cluster calculation manages to reproduce the q-dependence of the intensity well, but still 
fails to correctly account for the width of the experimentally observed feature at high q, 
thus signalling the impact of finite size effects, even in a 16 site cluster. 

The analysis of the momentum dependence of the loss function allows a more precise 
determination of the values of the model parameters. The best agreement is achieved with 
transfer integrals which mirror those derived from the band stucture calculation. In partic- 
ular, the inter-ladder hopping, txy is found to be 0.012 eV, thus confirming the quasi one- 
dimensional nature of the electronic system in a'-NaV205. Finally, the inter-site Coulomb 
interactions V^, Vb and Vxy are found to be of similar magnitude. These interactions then 
drive the electronic system close to a quantum critical point between a unordered state at 
room temperature and a zig-zag ordered state@ observed at low temperatures. 
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TABLES 



Ref. ta h tjyy U Va Vb Vxy 

Smolinski et al.i 
Horsch et al.il 
Damascelli et ahllZl 
Nishimoto et al.il 
Popovic et al.Ei 

Cuoco et al.il 0.4 0.2 0.15 4.0 0.8 0.8 0.9 

Sa et alS 0.38 0.17 0.012 2.8 0.37 0.10 0.43 

this work 0.38 0.17 0.012 2.8 0.8 0.6 0.9 

TABLE I. Literature values for the hopping parameters, t, the on-site Coulomb interaction U 
and the inter-site Coulomb interactions Va compared with the results of this work. 



ta 


tb 


txij 


U 


Va 


Vb 


0.38 


0.17 


0.012 


2.8 






0.35 


0.15 


0.3 


4.0 






0.3 


0.2 










0.3 


0.14 


0.05 


4.0 
6.82 




0.5 


0.4 


0.2 


0.15 


4.0 


0.8 


0.8 



E (eV) OxX Oyy 

1.1 V Sd-^xy V Sd^xy 

1.5 V 3d^xy V 3dixy 
3.1 2p^y3d^xy 

3.7 2p^Y 3dixy 02p^Y Sdi^y 

2-3 V Sd^xy * V 3dxz,yz V 3d^xy ^ V 3dxz,yz 



TABLE II. Energy positions and character of the electronic transitions contributing to the 
optical conductivities axx and ayy of Q;'-NaV205 in the low energy range. 
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FIGURES 



FIG. 1. Crystal structure of a'-NaV205. Vanadium atoms are shown as black spheres, oxygen 
atoms as grey spheres, (a) Layers of distorted VO5 pyramids (upper part) seperated by sodium 
ions (dark grey spheres) and the atomic ordering of a single vanadium oxide layer (lower part). 

(b) Projection of a single VO-layer onto the ab-plane, excluding the apex oxygens. The ladder 
structure is illustrated by the dashed lines. 

FIG. 2. Schematic structure of the ab planes in a'-NaV205 where the circles stand for the 
vanadium sites. The black sites define the two clusters used in the calculation. 

FIG. 3. Electron energy-loss spectra of a'-NaV205 measured with momentum transfer 
q = 0.1 A^^ parallel to the crystallographic a (solid line) and the crystallographic b direction 
(dashed line). 

FIG. 4. Spin projected densities of (a) the 2p states of the three inequivalent oxygen atoms 
and (b) the V 3d states in Q!'-NaV205 calculated for UeS=3 eV and an AFM order of V magnetic 
moments along the b direction. For comparison, the densities of the V 3d states calculated with a 
vector describing the spiral spin structure of Uy = 0.3 are also shown by dashed lines in (b). All 
energies are given relative to the valence band maximum. In the majority-spin V 3d DOS the two 
peaks just below the Fermi level and at 1 eV originate from a bonding and antibonding combina- 
tion of the 3dxy orbitals of the two V atoms on the same rung, respectively. The corresponding 
combinations of the minority-spin dxy states are shifted to 1.5 and 2.5 eV by the effective Coulomb 
repulsion. 

FIG. 5. The optical conductivity of a'-NaV205. (a) and (b) dyy from the EELS experiment. 

(c) and (d) show the same quantities from the LSDA and LSDA+?7 band-structure calculations with 
^^eff = (dotted line), 2 (solid line) and 3 eV (dashed line). The theoretical curves are broadened 
with a Lorentzian of width 0.2 eV. The vertical dotted lines illustrate the correspondence of the 
peak positions. 

FIG. 6. Decomposition of the optical conductivity of Q;'-NaV205 derived from the LSDA-|-i7 
calculations into contributions arising from transitions involving different initial and final states. 
In each case panels (a) and (c) show axx and (b) and (d), ayy. The thin solid lines indicate the 
total optical conductivities. In (a) and (b) the grey shading indicates initial states corresponding 
to occupied V 3dxy bands. In (c) and (d) light grey (dark grey) shading indicates antibonding 
V 3d'ixy (bonding V 3dixy) final states. The hatched area shows the contribution from interband 
transitions into higher lying empty V 3d states. 

FIG. 7. Optical conductivity of Q!'-NaV205 calculated for C/cff = 3 eV and different values of 
Uy, which defines the magnetic structure along the b direction (for details see text). 
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FIG. 8. Panels (a) and (c) show the loss function of a'-NaV205 measured using EELS with 
momentum transfer aligned along the a and b directions, respectively. Panels (b) and (c) show 
the corresponding calculated loss functions from a quarter filled t — J — V model, which have been 
broadened with a Gaussian function of width 0.3 eV. For the q = 0.1 calculated spectra, an 
unbroadened version has also been plotted to enable indentification of the individual plasmons. 
For the parameters used in the calculation, see text. 
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